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C"| , Abstract. We study the effect of tidal torques on the collapse of density peaks through 
^ ' the equations of motion of a shell of barionic matter falling into the central regions 
,H \ of a cluster of galaxies. We calculate the time of collapse of the perturbation taking 
d • into account the gravitational interaction of the quadrupole moment of the system with 



the tidal field of the matter of the neighbouring proto-clusters. We show that within 
high-density environments, such as rich clusters of galaxies, tidal torques slow down the 
collapse of low-u peaks producing an observable variation in the time of collapse of the 
shell and, as a consequence, a reduction in the mass bound to the collapsed perturbation. 
Moreover, the delay of the collapse produces a tendency for less dense regions to accrete 
less mass, with respect to a classical spherical model, inducing a biasing of over-dense 
regions toward higher mass. Finally we calculate the bias coefficient using a selection 
function properly defined showing that for a Standard Cold Dark Matter (SCDM) model 
this bias can account for a substantial part of the total bias required by observations on 
cluster scales. 
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1. Introduction 

It has long been speculated on the fundamental role that the angular momentum could 
play in determining the fate of collapsing proto-structures and several models have been 
proposed to correlate the galaxy type with the angular momentum per unit mass of the 
structure itself (Faber 1982; Kashlinsky 1982; Fall 1983). Some authors (see Barrow & 
Silk 1981; Szalay & Silk 1983 and Peebles 1990) have proposed that non-radial motions 
would be expected within a developing proto-cluster due to the tidal interaction of the 
irregular mass distribution around them, typical of hierarchical clustering models, with 
the neighbouring proto-clusters. The kinetic energy of these non-radial motions prevents 
the collapse of the proto-cluster, enabling the same to reach statistical equilibrium before 
the final collapse (the so-called previrialization conjecture by Davis & Peebles 1977, Pee- 
bles 1990). This effect may prevent the increase in the slope of the mass autocorrelation 
function at separations given by £(r, t) ~ 1, expected in the scaling solution for the rise 
in £(r, t) but not observed in the galaxy two-point correlation function. The role of non- 
radial motions has been pointed by several authors (see Davis & Peebles 1983: Gorski 
1988; Groth ct al. 1989; Mo et al. 1993; van de Weygaert & Babul 1994; Marzke et al. 
1995 and Antonuccio-Delogu & Colafrancesco 1995). Antonuccio-Delogu & Colafrancesco 
derived the conditional probability distribution f p k{v\i>) of the peculiar velocity around 
a peak of a Gaussian density field and used the moments of the velocity distribution to 
study the velocity dispersion around the peak. They showed that regions of the proto- 
clusters at radii, r, greater than the filtering length, Rf, contain predominantly non-radial 
motions. 

Non-radial motions change the energetics of the collapse model by introducing another 
potential energy term. In other words one expects that non-radial motions change the 
characteristics of the collapse and in particular the turn around epoch, t m , and conse- 
quently the critical threshold, 5 C , for collapse. Here, we want to remind that t m is the 
time at which the linear density fluctuations, that generate the cosmic structures, detach 
from the Hubble flow. The turn-around epoch is given by: 
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where pb is the mean background density, z is the redshift and 8 is the mean over- 
density within the non-linear region. After the turn around epoch, the fluctuations start 
to recollapse. As known for a spherical top hat model, the perturbation of the density 
field is completely collapsed when 

5 = 4 = (3/5)(^) 2 / 3 ^ 1.68 (2) 

where T c is the time of collapse which is twice the turn around epoch. One expects 
that non-radial motions produce firstly a change in the turn around epoch, secondly a new 
functional form for 5 C , thirdly a change in the mass function calculable with the Press- 
Schechter (1974) formula and finally a modification of the two-point correlation function. 
As we shall show in a forthcoming paper (Del Popolo & Gambera 1997b) non-radial mo- 
tions can reduce several discrepancies between the SCDM model and observations: the 
strong clustering of rich clusters of galaxies (£ C c(?~) — (r/25h~ 1 Mpc)~ 2 ) far in excess of 
CDM predictions (Bahcall & Soneira 1983), the X-ray temperature distribution function 
of clusters over-producing the observed cluster abundances (Bartlett & Silk 1993). 
For the sake of completeness, we remember that alternative models with more large-scale 
power than SCDM have been introduced in order to solve the latter problem. Several 
authors (Peebles 1984; Efstathiou et al. 1990; Turner 1991; White et al. 1993) have low- 
ered the matter density under the critical value (Cl m < 1) and have added a cosmological 
constant in order to retain a flat universe (Cl m + fl\ = 1) . The spectrum of the matter 
density is specified by the transfer function, but its shape is affected because of the fact 
that the epoch of matter-radiation equality is earlier, 1 + z eq being increased by a factor 
1/0 TO . Around the epoch z\ the growth of the density contrast slows down and ceases 
after z\. As a consequence the normalisation of the transfer function begins to fall, even 
if its shape is retained. 

Mixed dark matter models (MDM) (Bond et al. 1980; Shafi & Stecker 1984; Valdarnini 

6 Bonometto 1985; Holtzman 1989; Schaefer 1991; Schaefer & Shafi 1993; Holtzman & 
Primack 1993) increase the large-scale power because free-streaming neutrinos damp the 
power on small scales. Alternatively changing the primeval spectrum several problems of 
SCDM arc solved (Ccn et al. 1992). Finally, it is possible to assume that the threshold 
for galaxy formation is not spatially invariant but weakly modulated (2% — 3% on scales 
r > I0h~ 1 Mpc) by large scale density fluctuations, with the result that the clustering on 
large-scale is significantly increased (Bower et al. 1993). 

Moreover, this study of the role of non-radial motions in the collapse of density pertur- 
bations can help us to give a deeper insight in to the so-called problem of biasing. As 
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pointed out by Davis et al. (1985), unbiased CDM presents several problems: pairwise 
velocity dispersion larger than the observed one, galaxy correlation function steeper than 
that observed (see Liddle & Lyth 1993 and Strauss & Willick 1995). The remedy to these 
problems is the concept of biasing (Kaiser 1984), i.e. that galaxies are more strongly 
clustered than the mass distribution from which they originated. The physical origin of 
such biasing is not yet clear even if several mechanisms have been proposed (Rees 1985; 
Dekel & Silk 1986; Dekel & Rees 1987; Carlberg 1991; Cen & Ostriker 1992; Bower et 
al. 1993; Silk & Wyse 1993). Recently Colafrancesco, Antonuccio-Delogu & Del Popolo 
(1995, hereafter CAD) and Del Popolo & Gambera (1997a) have shown that dynamical 
friction delays the collapse of \ow-v peaks inducing a bias of dynamical nature. Because 
of the dynamical friction, under-dense regions in clusters (the clusters outskirts) accrete 
less mass than that accreted in absence of this dissipative effect and as a consequence 
over-dense regions are biased toward higher mass (Antonuccio-Delogu & Colafrancesco 
1994 and Del Popolo & Gambera, 1996). Non-radial motions acts in a similar way to 
dynamical friction: they delay the shell collapse consequently inducing a dynamical bias 
similar to that produced by dynamical friction. This dynamical bias can be evaluated 
defining a selection function similar to that given in CAD and using Bardeen, Bond, 
Szalay & Kaiser (1986, hereafter BBKS) prescriptions. 

The methods used in this paper are fundamentally some results of the statistics of Gaus- 
sian random fields, the biased galaxy formation theory and the spherical model for the 
collapse of density perturbations. In particular, we calculate the specific angular momen- 
tum acquired by protoclusters and the time of collapse of protoclusters using the Gaussian 
random fields theory and the spherical collapse model following Ryden's (1988a, hereafter 
R88a) approach. The selection function that we introduce is general and obtained by the 
only hypothesis of Gaussian density field. The approach and the final result is totally 
different from BBKS selection function and similar to that of Colafrancesco, Antonuccio 
& Del Popolo (1995). Only the biasing parameter is obtained from a BBKS approximated 
formula. This choice will be clarified in the following sections of the paper. 
The plan of the paper is the following: in Sect. 2 we obtain the total specific angular mo- 
mentum acquired during expansion by a proto-cluster. In Sect. 3 we use the calculated 
specific angular momentum to obtain the time of collapse of shells of matter around 
peaks of density having v c = 2,3,4 and we compare the results with Gunn & Gott's 
(1972, hereafter GG) spherical collapse model. In Sect. 4 we derive a selection function 
for the peaks giving rise to proto-structures while in Sect. 5 we calculate some values for 
the bias parameter, using the selection function derived, on three relevant filtering scales. 
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Finally in Sect. 6 we discuss the results obtained. 
2. Tidal torques 

The explanation of galaxies spins gain through tidal torques was pioneered by Hoyle 
(1949) in the context of a collapsing protogalaxy. Peebles (1969) considered the process 
in the context of an expanding world model showing that the angular momentum gained 
by the matter in a random comoving Eulerian sphere grows at the second order in pro- 
portion to i 5 / 3 (in a Einstein-de Sitter universe), since the proto-galaxy was still a small 
perturbation, while in the non-linear stage the growth rate of an oblate homogeneous 
spheroid decreases with time as t~ x . 

More recent analytic computations (White 1984; Hoffman 1986 and R88a) and numerical 
simulations (Barnes & Efstathiou 1987) have re-investigated the role of tidal torques in 
originating galaxies angular momentum. In particular White (1984) considered an anal- 
ysis by Doroshkevich (1970) that showing as the angular momentum of galaxies grows 
to first order in proportion to t and that Peebles's result is a consequence of the spher- 
ical symmetry imposed to the model. White showed that the angular momentum of a 
Lagrangian sphere does not grow either in the first or in the second order while the an- 
gular momentum of a non-spherical volume grows to the first order in agreement with 
Doroshkevich's result. 

Hoffman (1986) has been much more involved in the analysis of the correlation of the 
growth of angular momentum with the density perturbation 8{r). He found an angular 
momentum-density anticorrelation: high density peaks acquire less angular momentum 
than low density peaks. One way to study the variation of angular momentum with radius 
in a galaxy is that followed by R88a. In this approach the protogalaxy is divided into 
a series of mass shells and the torque on each mass shell is computed separately. The 
density profile of each proto-structure is approximated by the superposition of a spher- 
ical profile, 5(r), and a random CDM distribution, e(r), which provides the quadrupole 
moment of the protogalaxy. To the first order, the initial density can be represented by: 



being P{k) the power spectrum, while the density profile is (Ryden & Gunn 1987): 



p{v)= Pb [l + 5(r)\ [l + e{r)\ 



(3) 



where p b is the background density and e(r) is given by: 



(\e k \ 2 )=P(k) 



(4) 
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where v is the height of a density peak, £ (r) is the two-point correlation function, 7 
and i?* are two spectral parameters (BBKS, Eq. 4.6a, 4.6d) while d(^v, 7) is a function 
given in BBKS (Eq. 6.14). As shown by R88a the net rms torque on a mass shell centered 
on the origin of internal radius r and thickness Sr is given by: 

(Irl 2 ) 1 ^ = V30(yG) [{a 2m {rf){q 2m {rf) - (a 2m (r)q* 2m (r)) 2 ] 1/2 (6) 
where g; m , the multipole moments of the shell and a; m , the tidal moments, are given 

by: 

(g 2m (r) 2 > = -f-3M s \ / fc 2 dfcP (fc) J2 (fcr) 2 (7) 

(a 2m (r) 2 ) = /" dfcP (fc) j x (fcr) 2 (8) 

(a2m(r)^ m (r)) = ^PbM sh J kdkP(k)h (fcr) j 2 (fcr) (9) 

where M s ^ is the mass of the shell, ji(r) and j 2 {r) are the spherical Bessel function 
of first and second order while the power spectrum P{k) is given by: 

P(k) = Afc _1 [ln(l + 4.164fc)] 2 • (192.9 + 1340fc + 

+ 1.599 • 10 5 fc 2 + 1.78 • 10 5 fc 3 + 3.995 • 10 6 fc 4 )~ 1/2 (10) 

(Ryden & Gunn 1987). The normalization constant A can be obtained, as usual, by fixing 
that the mass variance at &h~ l Mpc, that is o" 8 , be equal to unity. Filtering the spectrum 
on cluster scales, Rf = 3h~ 1 Mpc 1 we have obtained the rms torque, r(r), on a mass shell 
using Eq. (g) then we obtained the total specific angular momentum, h(r,f), acquired 
during expansion integrating the torque over time (R88a Eq. 35): 

, T t~C 12 r (l-cosg) 3 Mpydfd 

h(r, V) — , / T-rr j- (11) 

mM sh h (tf-sintf) 4/3 /iW-/2(tf)| 
where the functions /i(t?), are given by R88a (Eq. 31) while the mean over- 

density inside the shell, S(r), is given by (R88a): 

3 r°° 

5( r ,u) - ^ / daa 2 6(a) (12) 
r Jo 

In Fig. 1 we show the variation of h(r, v) with the distance r for three values of the 
peak height v. The rms specific angular momentum, h(r,u), increases with distance r 
while peaks of greater v acquire less angular momentum via tidal torques. This is the 
angular momentum-density anticorrelation showed by Hoffman (1986). This effect arises 
because the angular momentum is proportional to the gain at turn around time, t m , 
which in turn is proportional to S(r, v)~i oc v~ z l 2 . 
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Fig. 1. The specific angular momentum, in units of Mq, Mpc and the Hubble time, to, for 
three values of the parameter v (y = 2 dotted line, v = 3 solid line, v — 4 dashed line) and for 
R f = 3/i _1 Mpc. 

3. Shell collapse time 

One of the consequences of the angular momentum acquisition by a mass shell of a 
proto-cluster is the delay of the collapse of the proto-structure. As shown by Barrow & 
Silk (1981) and Szalay & Silk (1983) the gravitational interaction of the irregular mass 
distribution of proto-cluster with the neighbouring proto-structures gives rise to non- 
radial motions, within the protocluster, which are expected to slow the rate of growth of 
the density contrast and to delay or suppress the collapse. According to Davis & Peebles 
(1977) the kinetic energy of the resulting non-radial motions at the epoch of maximum 
expansion increases so much to oppose the recollapse of the proto-structure. Numerical 
N-body simulations by Villumsen & Davis (1986) showed a tendency to reproduce this 
so-called previrialization effect. In a more recent paper by Peebles (1990) the slowdown 
in the growth of density fluctuations and the collapse suppression after the epoch of the 
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maximum expansion were re-obtained using a numerical action method. 
In the central regions of a density peak (r < 0.5i?/) the velocity dispersion attains nearly 
the same value (Antonuccio-Delogu & Colafrancesco 1995) while at larger radii (r > Rf) 
the radial component is lower than the tangential component. This means that motions 
in the outer regions are predominantly non-radial and in these regions the fate of the 
infalling material could be influenced by the amount of tangential velocity relative to the 
radial one. This can be shown writing the equation of motion of a spherically symmetric 
mass distribution with density n(r) (Peebles 1993): 

ft ft 71 (j 

_ nK) + + (2(^2) _ (u 2 }) _ + n[r) _ {vr) = (13) 

where (v r ) and (v^) are, respectively, the mean radial and tangential streaming veloc- 
ity. Eq. ( |l3| ) shows that high tangential velocity dispersion ((uf) > 2(v 2 )) may alter the 
infall pattern. The expected delay in the collapse of a perturbation may be calculated 
solving the equation for the radial acceleration (Peebles 1993): 

dVr = L 2 (r,u)_ 

dt Af 2 r 3 y{ ' { ' 

where L(r, v) is the angular momentum and g(r) the acceleration. Writing the proper 
radius of a shell in terms of the expansion parameter, a(ri, t): 

r(n,t) =na(n,t) (15) 

remembering that 
An 

M = —p b (n,t)a 3 (n,t)rf (16) 

and that pb = where Hq is the Hubble constant and assuming that no shell 

crossing occurs so that the total mass inside each shell remains constant, that is: 

/ ,n Pi(n,t) , . 

p{r - t) = ^-T) (17) 



Eq. (|14| ) may be written as: 
2 a _ H 2 (l + 5) ±G 2 L 2 



dt 2 2a 2 Hi(l + 6) 2 rj°a 3 

or integrating the equation once more: 



(18) 



2 I 



8G " L2 2 l-da-2C (19) 



where C is the binding energy of the shell. The value of C can be obtained using the 
condition for turn around 4j = when a = a max leading to the new equation: 
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~dt . 



Hi 



8G 2 L 2 



da 



(20) 



Equation ( |l4| ) or equivalently Eq. ([l^) may be solved using the initial conditions: 
(<M ) = 0> a = a max — 1/5 and using the function h(r } v) = L{r, u)/M s h found in Sect. 2 
to obtain the time of collapse, T c (r, v). 

In Figs. 2 ™ 4 we compare the results for the time of collapse, T c , for v = 2,3,4 with the 
time of collapse of the classical GG spherical model: 



T co{r, v) = — \5{r,v) 



-3/2 



(21) 




Fig. 2. The time of collapse of a shell of matter in units of the age of the Universe to for v = 2 
(dotted line) compared with GG's model (solid line). 

As shown the presence of non-radial motions produces an increase in the time of col- 
lapse of a spherical shell. The collapse delay is larger for a low value of v and becomes 
negligible for v > 3. This result is in agreement with the angular momentum-density anti- 
correlation effect: density peaks having low value of v acquire a larger angular momentum 
than high v peaks and consequently the collapse is more delayed with respect to high v 
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Fig. 3. The time of collapse of a shell of matter in units of the age of the Universe to for v = 3 
(dotted line) compared with GG's model (solid line). 

peaks. Given T c (r, v) we also calculated the total mass gravitationally bound to the final 
non-linear configuration. There are at least two criteria to establish the bound region to 
a perturbation S(r): a statistical one (Ryden 1988b), and a dynamical one (Hoffman & 
Shaham 1985). The dynamical criterion, that we have used, assumes that the binding 
radius is given by the condition that a mass shell collapse in a time, T c , smaller than the 
age of the universe to: 

T c (r,v)<t (22) 

We calculated the time of collapse of GG spherical model, T c o(r, v) : using the density 
profiles given in Eq. (|5|) for 1.7 < v < 4 and then we repeated the calculation taking into 
account non-radial motions obtaining T c (r,v). Then we calculated the binding radius, 
rb(y), for a GG model solving T c o(r, v) < to for r and for several values of u, while we 
calculated the binding radius of the model that takes into account non-radial motions, 
rb(j/), repeating the calculation, this time with T c (r, v) < to. We found a relation between 
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Fig. 4. The time of collapse of a shell of matter in units of the age of the Universe to for v = 4 
(dotted line) compared with GG's model (solid line). 

v and the mass of the cluster using the equation: M = ^r^pt,. 

In Fig. 5 we compare the peak mass obtained from GG model, using Hoffman & Sha- 
ham's (1985) criterion, with that obtained from the model taking into account non-radial 
motions. As shown for high values of v [y > 3), the two models give the same result for 
the mass while for v < 3 the effect of non-radial motions produces less bound mass with 
respect to GG model. 

Before concluding this section we want to discuss the applicability of the spherical model 
and consequently the use of spherical shells in our model. The question of the applicablity 
of the idealized spherical collapse model and that of the secondary infall model (SIM) to 
realistic systems and initial conditions is almost as old as the model itself (GG; Gunn 
1977; Filmore & Goldreich 1984; Quinn ct al. 1986; Zurek et al. 1988; Quinn & Zurck 
1988; Warren et al. 1991; Crone et al. 1994). In a recent paper Zaroubi et al. (1996) in- 
vestigate the applicability of the quoted model to realistic systems and initial conditions 
by comparing the results obtained using the spherical model and SIM with that of a set 
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Fig. 5. The mass accreted by a collapsed perturbation, in units of 1O 15 M , taking into account 
non-radial motions effect (dotted line) compared to GG's mass (solid line). 

of simulations performed with different N-body codes (Treecode and a " monopole term" 
code). They introduced a numerical model to trace the evolution of density peaks under 
the assumption of asphcrically symmetric force and realistic initial conditions instead 
of the spherically symmetric force and initial conditions as assumed in the SIM. They 
obtained a good agreement between the SIM and the simulations. 

4. Tidal field and the selection function 

According to biased galaxy formation theory the sites of formation of structures of mass 
~ M must be identified with the maxima of the density peak smoothed over a scale 
Rf (M oc R 3 f)- A necessary condition for a perturbation to form a structure is that it 
goes non-linearly and that the linearly extrapolated density contrast reaches the value 
5(r) > S c = 1.68 or equivalently that the threshold criterion v t > 5 c /ao(Rf) is satis- 
fied, being a (Rf) the variance of the density field smoothed on scale Rf. When these 




13 



conditions are satisfied the matter in a shell around a peak falls in toward the cluster 
center and virializes. In this scenario only rare high v peaks form bright objects while 
low v peaks (y ~ 1) form under-luminous objects. The kind of objects that forms from 
non-linear structures depends on the details of the collapse. Moreover, if structures form 
only at peaks in the mass distribution they will be more strongly clustered than the mass. 
Several feedback mechanisms have been proposed to explain this segregation effect (Rees 
1985; Dekel & Rees 1987). Even if these feedback mechanisms work one cannot expect 
they have effect instantaneously, so the threshold for structure formation cannot be sharp 
(BBKS). To take into account this effect BBKS introduced a threshold or selection func- 
tion, t(v/v t ). The selection function, t{v/vt), gives the probability that a density peak 
forms an object, while the threshold level, v tl is defined so that the probability that a 
peak forms an object is 1/2 when v = v t . The selection function introduced by BBKS 
(Eq. 4.13), is an empirical one and depends on two parameters: the threshold v t and the 
shape parameter q: 

" <») 

If q — ► oo this selection function is a Heaviside function d(v — v t ) so that peaks with 
v > v t have a probability equal to 100% to form objects while peaks with v < u t do not 
form objects. If q has a finite value sub-^t peaks are selected with non-zero probability. 
Using the given selection function the cumulative number density of peaks higher than v 
is given, according to BBKS, by: 

/oo 
t(v/v t )N pk (v)dv (24) 

where N p k (v) is the comoving peak density (see BBKS Eq. 4.3). A form of the selection 
function, physically founded, can be obtained following the argument given in CAD. In 
this last paper the selection function is defined as: 



p OO 

v) = J p\5,(S(rMt,v)),(?-s(rMt,v)]d6 (25) 



where the function 



P p, (J(r))] ^ exp (-™P) (26) 

gives the probability that the peak overdensity is different from the average, in a 
Gaussian density field. The selection function depends on v through the dependence of 
S(r) on v. As displayed, the integrand is evaluated at a radius rut which is the typical 
radius of the object that we are selecting. Moreover, the selection function t(y) depends 
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on the critical overdensity threshold for the collapse, S c , which is not constant as in a 
spherical model (due to the presence, in our analysis, of non-radial motions that delay 
the collapse of the proto-cluster) but it depends on v. The dependence of S c on v can be 
obtained in several ways, for example according to Peebles (1980) the value of 8 C depends 
on the ratio T c /t m between the perturbation collapse time, T c , and its turn around time, 

<-mr 

Non-radial motions slow down the collapse of the mass shell with respect to the 
GG collapse time changing the value of T c . Using the calculated time of collapse for a 
given shell, and its dependence on i>, 8 c {y) can be calculated using Eq. (^). An analityc 
determination of S c (i>) can be obtained following a technique similar to that used by 
Bartlett & Silk (1993). Using Eq. (|l9| ) it is possible to obtain the value of the expansion 
parameter of the turn around epoch, a max , which is characterized by the condition 4| = 0. 
Using the relation between v and Si, in linear theory (Peebles 1980), we can find C that 
substituted in Eq. ( |l9| ) gives at turn around: 

8G 2 f L 2 da] 



8 c (v) = 5 c0 



1 



(28) 



n 3 o Hgr} 5(l + S) 2 

where 5 c q — 1.68 is the critical threshold for GG's model. In Fig. 6 we show the 
overdensity threshold as a function of v. As shown, 5 c (v) decreases with increasing v. 
when v > 3 the threshold assumes the typical value of the spherical model. This means, 
according to the cooperative galaxy formation theory, (Bower et al. 1993) that structures 
form more easily if there are other structures nearby, i.e. the threshold level is a decreasing 
function of the mean mass density. Known 5 c {v) and chosen a spectrum, the selection 
function is immediately obtainable through Eq. (|2^) and Eq. (p(f). The result of the 
calculation, plotted in Fig. 7, for two values of the filtering radius, {Rf = 3, 4 h^Mpc), 
shows that the selection function, as expected, differs from an Heaviside function (sharp 
threshold). The shape of the selection function depends on the values of the filtering 
length Rf and on non-radial motions. The value of v at which the selection function 
t(v) reaches the value 1 {t(v) ~ 1) increases for growing values of the filtering radius, 
Rf. This is due to the smoothing effect of the filtering process. The effect of non-radial 
motions is, firstly, that of shifting t(i>) towards higher values of v, and, secondly, that 
of making it steeper. The selection function is also different from that used by BBKS 
(Tab. 3a). Finally it is interesting to note that the selection function defined by Eq. ( |25| ) 
and Eq. ( p6~l) is totally general, it does not depend on the presence or absence of non- 
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radial motions. The latter influences the selection function form through the change of 
5 C induced by non-radial motions itself. 

5. The bias coefficient 

A model of the Universe in which light traces the mass distribution accurately (unbiased 
model) is subject to several problems. As pointed out by Davis et al. (1985) an unbiased 
CDM produces a galaxy correlation function which is steeper than that observed and a 
pairwise velocity dispersion larger than that deduced from redshift surveys. A remedy to 
this problem can be found if we do not assume that light traces mass and adopt the biasing 
concept, i.e., galaxies are more clustered than the distribution of matter in agreement 
with the concept of biasing inspired by Kaiser (1984). The observations show that clusters 
of galaxies cluster more strongly than galaxies, in the sense that the dimensionless two- 
point correlation function, £ cc (r), is much larger than the galaxy two-point function, 
£ gg (r). The galaxy two-point correlation function £ gg {r) is a power-law: 

«» w =fey <29) 
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Fig. 6. The selection function, t(y), for Rf — 3h~ Mpc (the solid line plots the selection 
function obtained without taking into account the effects of non-radial motions; the dotted 
line plots the selection function obtained taking into account the effects of non-radial motions) 
and for 4/i -1 Mpc (the short dashed line plots the selection function obtained without taking 
into account the effects of non-radial motions; the long dashed line plots the selection function 
obtained taking into account the effects of non-radial motions) . 

with a correlation length ro. 3 ~ 5/i _1 Mpc and a slope 7 ~ 1.8 for r < lOh^ 1 Mpc 
(Davis & Peebles 1983; Davis et al. 1985; Shanks et al. 1989), (some authors disagree with 
this values; for example Strauss et al. 1992 and Fisher et al. 1993 find ro, 9 ~ 3.79/i _1 Mpc 
and 7 ~ 1.57). As regards the clusters of galaxies the form of the two-point correlation 
function, £ cc (r), is equal to that given by Eq. (^9|). Only the correlation length is different. 
In the case of clusters of galaxies the value of ro jC is uncertain (see Bahcall & Soneira 
1983; Postman et al. 1986; Sutherland 1988; Bahcall 1988; Dekel et al. 1989; Olivier et 
al. 1990 and Sutherland & Efstathiou 1991) however it lies in the range ro jC — 12~25h^ 1 
Mpc in any case larger than ro, g . One way of defining the bias coefficient of a class of 
objects is that given by (BBKS): 
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b(R f ) = ?L + 1 (30) 



where (£) is: 

7 6» 



l- 7 2 



/ ( ~ ) a;,,, .(/'),//' (.tii 



We want to remember that, as shown by Coles (1993), the biasing parameter can be 
also estimated by means of the ratio of the amplitudes of the correlation function, £(r), 
and the matter auto-covariance function, T(r): 

b 2 (r) = ^ (32) 

or by means of the ratio of the cumulative integral of the two-points correlation 
function (K^(r) = £(r )r 2 dr) and that of the auto-covariance function (Ja(r) = 
f,nr)r 2 dr): 

Mr) = Mr) (33) 
or finally the ratio of galaxy, Q(k), to mass, P(k), power spectra: 

*M - (34) 

As stressed by Coles (1993) a local bias generally produces a different response in each 
of these descriptors. Even the qualitative behaviour of the limit of large scales can be 
different, i.e. b(r) can increase or decrease. So one should decide very carefully which one 
of these definitions must be used, when discussing the behaviour of the biasing parameter. 
Being conscious of these difficulties we have chosen one of the most popular descriptor of 
the biasing parameter (see also Lilje 1990, Liddle & Lith 1993, Croft & Efstathiou 1994, 
CAD) in order to make comparisons with other models. From Eq. ( |3l| ) it is clear that 
the bias parameter can be calculated once a spectrum, P(k), is fixed. The bias parameter 
depends on the shape and normalization of the power spectrum. A larger value is obtained 
for spectra with more power on large scale (Kauffmann et al. 1996). In this calculation we 
continue to use the standard CDM spectrum (Oo = 1, h = 0.5) normalized imposing that 
the rms density fluctuations in a sphere of radius &h~ l Mpc is the same as that observed 
in galaxy counts, i.e. as = a(8h~ 1 AIpc) = 1. The calculations have been performed for 
three different values of the filtering radius (Rf — 2, 3, 4 h~ 1 Mpc). The values of b, that 
we have obtained, are respectively, in increasing order Rf, 1.6, 1.93 and 2.25. 
As shown, the value of the bias parameter tends to increase with Rf due the filter effect 
of t(i>). As shown t{v) acts as a filter, increasing the filtering radius, Rf, the value of v 



18 



at which t(u) ~ 1 increases. In other words when Rf increases, t(u) selects density peaks 
of a larger height. The reason for this behavior must be searched in the smoothing effect 
that the increasing of the filtering radius produces on density peaks. When Rf increases 
the density field smooths and t{v) has to shift towards a higher value of v in order to 
select a class of object of fixed mass M. 

6. Conclusions 

In this paper we have studied the role of non-radial motions on the collapse of density 
peaks solving numerically the equations of motion of a shell of barionic matter falling 
into the central regions of a cluster of galaxies. We have shown that non-radial motions 
produce a delay in the collapse of density peaks having a low value of v while the collapse 
of density peaks having v > 3 is not influenced. A first consequence of this effect is a 
reduction of the mass bound to collapsed perturbations and an increase in the critical 
threshold, 5 C , which now is larger than that of the top-hat spherical model and depends 
on v. This means that shells of matter of low density have to be subject to a larger 
gravitational potential, with respect to the homogeneous GG's model, in order to collapse. 
The delay in the proto-structures collapse gives rise to a dynamical bias similar to that 
described in CAD whose bias parameter may be obtained once a proper selection function 
is defined. The selection function found is not a pure Hcavisidc function and is different 
from that used by BBKS to study the statistical properties of clusters of galaxies. Its 
shape depends on the effect of non-radial motions through its dependence on 8 c {v). The 
function t{v) selects higher and higher density peaks with increasing value of Rf due to 
the smoothing effect of the density field produced by the filtering procedure. Using this 
selection function and BBKS prescriptions we have calculated the coefficient of bias b. 
On clusters scales for Rf = Ah~ 1 Mpc we found a value of b = 2.25 comparable with 
that obtained from the mean mass-to-light ratio of clusters, APM survey, or from N- 
body simulations combined with hydrodynamical models (Frenk et al. 1990). Morever, 
the value of the coefficient of biasing b that we have calculated is comparable with the 
values of b given by Kauffmann et al. (1996). This means that non-radial motions play a 
significant role in determining the bias level. In our next paper (Del Popolo & Gambera, 
in preparation) we make a more detailed analysis on the problem of the bias. 
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